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ABSTRACT 

Optimum frequency for propagation of sound in shallow 
sound channels was studied using two acoustic transmission 
loss models. The split-step Parabolic Equation (PE) model, 
a full-wave model; and the Fast Asymptotic Coherent 
Transmission loss model, version 9H (FACT 9H) , a ray- tracing 
model; were tested against experimental data collected by 
Dosso and Chapman (1984) in the northeast Pacific Ocean. 
The models were found to be valid predictors of optimum 
frequency for the shallow sound channel observed by Dosso 
and Chapman. Both models were then used to predict optimum 
frequency for two sound .velocity profiles obtained in a 
high- latitude deep ocean basin under summer conditions, 
exhibiting shallow sound channels. As expected, the split- 
step PE model adequately predicted optimum frequencies for 
these cases. The FACT 9H model did not produce reasonable 
results for optimum frequencies. 
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PREFACE 



The Environmental Acoustic Research Group at the Naval 
Postgraduate School(NPS) is engaged in research to establish 
beneficial and detrimental environmental effects important 
to present and future Navy acoustics systems. 

Pursuant to the above objectives, environmental and 
acoustics models are used to interpret and predict the 
complex results obtained when actual experimental or opera- 
tional scenarios are utilized. The present work is part, of 
a series of investigations concerning shallow sound channels 
in ASW [Ref. 1] . 

The Environmental Acoustic Research Group presently 
consists of 7 professors from the Departments of 
Oceanography, Physics and Electrical and Computer 
Engineering. About 10 NPS graduate students are partici- 
pating in the research group. 
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I. INTRODUCTION 



A. BACKGROUND 

Despite technological advances in ship silencing, the 
acoustic signature produced by a submarine's main propulsion 
system and auxiliary machinery remains the single most 
easily exploited tool for detecting, tracking and localizing 
the submarine. The transmission of sound through the ocean 
is predictable given a sufficiently well-defined knowledge 
of the environment, the source and the receiver. 
Transmission phenomena can become extremely complex, even 
for relatively simple sound velocity profiles. 

Acoustic transmission ' loss models have been developed 
for implementation on digital computers. The models gener- 
ally take one of two approaches: either solving the 
acoustic wave equation or employing ray-tracing techniques. 
Both approaches require assumptions, simplifications and 
approximations. For real situations some models will 
predict transmission loss more accurately than others 
depending on the computational technique involved, the 
method of solution and the extent to which the model accepts 
real conditions for boundaries at the sea surface, the sea 
floor and between different water masses. 

The situation of interest here is one of the more 
complicated cases: the shallow sound channel. This situ- 
ation develops in higher latitudes of the world's oceans, 
where precipitation is heavy and evaporation is small. A 
sound velocity profile reveals the presence of a mixed layer 
resulting from winds and turbulence which may extend from 
the surface to depths of several tens of meters. Below the 
mixed layer a shallow sound channel (SSC) is formed by a 
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combination of temperature and salinity gradients, resulting 
in a sound speed minimum (channel axis) at roughly 80 to 200 
meters. The lower edge of the SSC is bounded by an area of 
rapidly increasing sound speed. Below the SSC the combined 
effects of decreasing temperature and increasing pressure 
can also result in a second or deep sound channel (DSC) with 
an axis at about 800 to 1000 meters in mid- latitudes . Sound 
channels act as acoustic waveguides, trapping certain 
frequencies and allowing them to propagate especially well, 
while other frequencies "leak" out of the channel and propa- 
gate poorly. A more complete discussion of the theoretical 
background for this paper is contained in Chapter II. 

The various layers and channels discussed above occupy 
the top several hundred to one thousand meters of the 
oceans. Submarines operate within these depths. These 
features of the sound velocity profile strongly affect the 
propagation of the submarine's acoustic signature, and the 
active and passive operation of the submarine's sonar equip- 
ment and that of its hunters. Antisubmarine warfare (ASW) 
operators must understand these effects and have some means 
of predicting them. With reliable predictions operators 
will be able to select frequencies and modes of operation to 
maximize the effectiveness of their sonar equipment in 
detecting and tracking the submarine. 

B . OBJECTIVE 

This study seeks to use two acoustic transmission loss 
models to predict the behavior of sound in typical shallow 
sound channels. Both models are available at the Naval 
Postgraduate School (NPS), accessed through the IBM 3033 
computer . 

The first model, the split-step Parabolic Equation (PE) 
model, uses a parabolic approximation to the basic wave 
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equation and solves for pressure levels as a function of 
range. It is considered to be a low-frequency (below 800 
Hz) model, since higher frequencies require large amounts of 
computational time. It is not commonly used in support of 
fleet operations due to the large computational times 
required . 

The second model, the Fast Asymptotic Coherent 
Transmission loss model version 9H (FACT 9H) , employs ray- 
tracing techniques to predict the amount of acoustic energy 
arriving at a receiver as a function of range. This model 
is less rigorous in predicting the effects of different 
types of boundaries and it does not allow for variations of 
the sound velocity profile with range. It is simpler and 
requires less computational time. FACT is generally consid- 
ered more appropriate than the PE model for frequencies 
above 800 Hz. The FACT model is widely used in acoustic 
transmission loss predictions provided to fleet users in the 
U. S. Navy. Both the PE model and FACT 9H are available at 
the Naval Postgraduate School, accessed through the IBM 3033 
computer . 

In 1984 Dosso and Chapman [Ref. 2] reported an experi- 
ment on propagation loss that they had conducted in the 
northeast Pacific Ocean. Sound velocity profiles indicated 
the presence of a shallow sound channel overlying a deep 
sound channel. They placed sources in the shallow sound 
channel, with receivers in the shallow sound channel and the 
deep sound channel. They measured propagation losses over 
ranges extending from zero to 35 km for propagation both 
within the SSC and across the boundary with the DSC. Using 
a split-step parabolic equation model, Dosso and Chapman 
were able to obtain good agreement between their experi- 
mental results and the predictions of the model. A full 
discussion of Dosso and Chapman's experiment is contained in 
Chapter III. 
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This study examines the Dosso and Chapman experiment by 
using the Naval Postgraduate School’s PE and FACT 9H models. 
Chapter IV compares these predictions with the Dosso and 
Chapman's experimental results, supporting the use of the 
two models to predict the optimum frequency for sound propa- 
gation within a shallow sound channel. The models are then 
applied to two sound velocity profiles of interest to the 
Naval Postgraduate School's Environmental Acoustic Research 
Group. These profiles were obtained under summer conditions 
in a deep ocean basin (4000 meters) at high latitude. 
Chapter V describes the results of the NPS PE model as 
applied to these profiles. Chapter VI discusses the FACT 9H 
model and its predictions for the Dosso and Chapman profiles 
and the high- latitude summer profiles. Chapter VII then 
summarizes this study's results and conclusions and presents 
recommendations for further investigation. 
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II. THEORY 



A. SOUND IN THE SEA 

The manner in which sound is transmitted through a 
medium depends on the velocity with which it propagates. 
The velocity of sound at any given point in the ocean is 
calculated from the temperature, salinity and water depth. 
Since pressure varies linearly with water depth, these two 
quantities can be treated as interchangeable. In shallow 
depth temperature is generally the dominant factors 
affecting the speed of sound, but in the deep depth pressure 
is. Thus the velocity of sound in the ocean varies in all 
dimensions: horizontally, vertically and temporally 
[Ref . 3] . 

The ocean may be thought of as composed of many layers 
of water, in each of which sound has a different velocity. 
Where the sound velocity profile (SVP) shows a velocity 
minimum, sound is refracted toward that depth, resulting in 
a focusing of energy at that velocity minimum. This depth 
is often referred to as a sound channel axis. 

A typical sound velocity profile like that of Figure 2.1 
[Ref. 4] shows a surface layer of seawater which is subject 
to the actions of sun, wind and precipitation. This layer 
is often well-mixed, but its temperature, salinity and 
thickness may vary greatly from place to place and from day 
to day. Typically the surface layer shows a constant or 
slightly increasing sound velocity profile. The surface 
layer may extend to a depth of several tens of meters. At 
some point below the surface layer a strong thermal gradient 
may show up. Temperature and therefore sound velocity 
decrease sharply with depth, creating a seasonal thermo- 
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cline. Deeper still lies the main thermocline, formed by 
another strong negative thermal gradient. Below this is the 
deep isothermal layer, which extends to the sea floor. In 
this layer salinity is essentially constant, as is tempera- 
ture. Pressure and velocity increase very nearly linearly 
with depth. The layers bounded above by the main thermo- 
cline and below by the increasing velocity gradient are 
called the deep sound channel (DSC). Sound can be trapped 
in this channel and refracted repeatedly toward the velocity 
minimum, which is called the DSC axis. This axis may occur 
at depths ranging from several hundred to over 1000 meters. 

Certain environmental conditions can cause a sound 
velocity profile like that shown in Figure 2.2. A combina- 
tion of winds, high precipitation and low evaporation such 
as is found in the northeastern Pacific Ocean can lead to 
the formation of a secondary or shallow sound channel (SSC) 
just below the seasonal thermocline. Just as sound can be 
trapped in the DSC sound can also be trapped in the SSC 
[Ref. 2]. 

B . ACOUSTIC WAVEGUIDES 

Sound channels act as acoustic waveguides. The ray 
theory of sound propagation predicts that sound can be 
refracted and trapped in ducts. Trapping depends on the 
thickness of the duct and the wavelength of the sound; if 
the wavelength is "too large" to be contained within the 
vertical dimension of the duct, the sound will propagate 
poorly. Frequencies corresponding to those longer wavel- 
engths will propagate poorly, if at all; higher frequencies 
with shorter wavelengths will propagate much better. Normal 
mode theory provides a tool for understanding these propaga- 
tion effects [Ref. 4]. The transition from poor to good 
propagation is marked by the cutoff frequency. This 
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frequency does not mark an abrupt end to propagation below 
which sound cannot be transmitted at all; rather it is a 
lower boundary below which propagation is less effective. 

Acoustic propagation in a sound channel is exactly anal- 
ogous to the propagation of radio waves in a ground-based 
duct. Urick [Ref. 3] provides an equation from the theory 
of radio propagation for the maximum wavelength which will 
"fit" in such a duct. Dosso and Chapman [Ref. 2] use this 

equation to calculate the cutoff frequency for their 

80-meter thick shallow sound channel as about 250 Hz. 

Urick gives another equation derived from radio propaga- 
tion theory to calculate the cutoff frequency: [Ref. 3] 

| 

X v != 4.7 x 10" 3 * H 3/2 (2.1) 

IHglX 

where ^ax the max i mum wavelength of the trapped sound 

in feet, and H is the thickness of the duct in feet. 
Converted : ato metric units and using a sound velocity of 
1483 m/s, jhis formula yields a cutoff frequency of about 
243 Hz for Dosso and Chapman's 80-meter channel. A third 

formula is given by Kinsler, Frey, Coppens and Sanders 
[Ref. 4] 

f co * 2 x 10 5 / D 3/2 (2.2) 

for D in meters. This yields a cutoff frequency of about 
280 Hz for this channel. All of these calculations agree 
within useful ranges, since the cutoff frequency does not 
mark an abrupt end to propagation. 

C. TRANSMISSION LOSS 

Several mechanisms act to decrease the energy of an 
acoustic wave as the wavefront moves farther away from its 
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source. The most obvious loss mechanism is spreading. Once 
an outwardly propagating wavefront reaches a reasonable 
distance from a sufficiently small source, the wavefront can 
be considered for all practical purposes to be spherical. 
Thus , the energy in the wave is spread out over an area that 
increases as the square of the range from the source. In 
decibels, spherical spreading alone would produce a trans- 
mission loss (TL) of 

TL , =20 log, « R (2.3) 
spher ^10 v ' 

If the sound is somehow contained so that spreading may 
only take place horizontally, as is the case when the sound 
source is inside a channel or duct, the spreading is cylin- 
drical rather than spherical and the loss is given by 

TL cyl = 10 log 10 R (2 ' A) 

Other mechanisms cause additional energy to be lost. 
These attenuation mechanisms include scattering and 
absorption. The degree of attenuation is strongly dependent 
on the frequency of the sound, with higher frequencies 
experiencing greater attenuation. Thorp [Ref. 5] developed 
extensive tables to describe the attenuation of sound as a 
function of frequency. 

Sound can be scattered from rough boundaries at the sea 
surface or bottom, or from discontinuities at the edges of 
the layers. Bubbles, marine organisms and suspended partic- 
ulate matter will also scatter sound [Ref. 4] . Scattering 
is generally more important for higher frequencies than for 
lower ones . 

Sound energy can also be lost through absorption. 
Viscosity effects in the interaction of water molecules as 
the sound wave passes account for one absorption mechanism. 



19 



Chemical reactions are also important sources of absorptive 
loss. Of the many chemical constituents of seawater, 
certain compounds appear to undergo continuous dissociation 
into ions and reassociation. This relaxation mechanism 
takes up acoustic energy. The most important absorption 
mechanism for frequencies near 1000 Hz appears to be the 
borate-boric acid relaxation mechanism [Ref. 6]. 

In general, attenuation increases logarithmically with 
frequency. The attenuation losses due to such mechanisms 
are so small at low frequencies that the measurements must 
be made over extremely long ranges to detect any losses. 
The precision with which the measurements must be made and 
the inhomogeneities in the ocean further complicate the 
task. Urick [Ref. 3] provides an excellent discussion of 
absorption and attenuation^ 

Two factors are operating with opposing effects on sound 
being propagated through sound channels. First, attenuation 
mechanisms tend to cause losses of acoustic energy. Second, 
the focusing effect of the channel tends to reduce the 
losses. In the ideal case, a sound channel bounded above 
and below by smooth interfaces between layers with no 
absorption or leakage should produce losses at about the 
rate due to cylindrical spreading alone. The minimum 
frequency which will be trapped well in a channel is the 
cutoff frequency. At frequencies near 1000 Hz, absorption 
mechanisms become important, and scattering losses also come 
into play. Therefore it should be possible to find some 
optimum frequency for transmission of sound in a channel of 
given thickness, lying somewhere above the cutoff frequency 
and below the frequency range where absorption and scat- 
tering losses become large. 
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Figure 2.1 Typical mid-latitude sound velocity profile. 
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Figure 2.2 Sound velocity profile 
showing shallow sound channel. 
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III. THE DOSSO AND CHAPMAN EXPERIMENT 



A. DATA COLLECTION 

Dosso and Chapman [Ref. 2] conducted an experiment in 
sound propagation within and below the shallow sound 
channel. During summer conditions in the northeast Pacific 
Ocean off the west coast of Canada, CFAV ENDEAVOUR took up a 
station and monitored data collection equipment. A single 
calibrated hydrophone was suspended below the vessel at a 
depth of 101 meters, approximately the shallow sound channel 
axis. An array of hydrophones was suspended at 417 meters 
near the deep sound channel axis. A second vessel, CVAF 
PARIZEAU, opened range from CFAV ENDEAVOUR, dropping explo- 
sive charges every 1.8 km (approximately one nautical mile). 
The charges were set to detonate at 98 meters, near the 
shallow sound channel axis. 

CFAV PARIZEAU also measured sound velocity profiles at 
16-kilometer intervals over the track (Figures 3.1 through 
3.3). Later analysis revealed that the sound channel deep- 
ened over the portion of the track between eight and twenty 
kilometers, so that the sources were not located within the 
SSC and sound trapping was greatly reduced for sources deto- 
nated between those ranges. Section D below discusses steps 
that Dosso and Chapman took to account for this deepening of 
the SSC. 

The received acoustic energy was measured as 1/3 octave 
band averages for each shot for both the shallow and deep 
hydrophones. Dosso and Chapman were able to distinguish 
bottom- interacting propagation paths by the different 
arrival times for shots in the first 35 kilometers of the 
track, and these arrivals were subtracted from the resulting 
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energy levels. Bottom interactions could not be isolated 
from more distant signals. Propagation loss was calculated 
as known source level minus the received level measured at 
the hydrophones . 

B. RESULTS AND DISCUSSION 

For the shallow sound channel [Ref. 2], measured propa- 
gation loss was plotted as a function of range for the 1/3 
octave bands. The propagation loss demonstrated strong 
dependence on range, as expected. An optimum frequency for 
propagation of about 800 Hz was observed. 

Figure 3.4 shows the propagation loss measured by Dosso 
and Chapman for both source and receiver located within the 
shallow sound channel. The solid curve represents geome- 
trical spreading loss, calculated from 

H = 10 log.. (R • R ) (3.1) 

1U o 

where range R and ocean depth R 0 are measured in meters. 
Figure 3.4 (a) represents propagation loss below the optimum 
frequency of 800 Hz; Figure 3.4 (b) represents propagation 
loss above 800 Hz. The close agreement between the geome- 
trical loss curve and the 800 Hz data over much of the range 
is noteworthy. 

Figure 3.5 displays the propagation loss within the 
shallow sound channel for the optimum propagation frequency 
of 800 Hz, with the effects of chemical absorption removed. 
Again the solid line represents geometrical spreading for 
comparison. Note the abrupt increase in losses over the 
8-20 km range. Over this range the sound channel axis deep- 
ened, so that the sources were no longer within the duct. 

Figure 3.6 presents the same data plotted as a function 
of frequency for three-point range - averages centered at 24 
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and 33 km. The solid circles represent actual measurements; 
the open circles represent values adjusted upward to remove 
the effects of chemical absorption. The dotted lines repre- 
sent the expected loss from geometrical spreading alone at 
800 Hz. This display shows more clearly the optimum 
frequency of 800 Hz and how closely the losses for this 
frequency approximate the losses due to geometrical 
spreading alone. It was clear that frequencies above and 
below 800 Hz experienced greater losses. Diffraction or 
"leakage” from the shallow sound channel accounted for the 
losses at lower frequencies. High frequencies experienced 
losses due to chemical absorption mechanisms and scattering. 

C. COMPARISON BETWEEN SSC AND DSC PROPAGATION 

The deep sound channel' trapped lower frequencies (longer 
wavelengths) than the shallow sound channel because the DSC 
is much thicker. Figure 3.7 demonstrates the effects of 
receiver depth on propagation loss. The sources were 

located at 98 meters, near the shallow sound channel axis. 
One receiver was also near the shallow sound channel axis at 
101 meters. The other receiver was positioned at 417 
meters, near the deep sound channel axis. At 160 Hz sound 
leaked easily out of the shallow sound channel but was well 
trapped within the deep sound channel. At 315 Hz, just 

above the cutoff frequency for the shallow sound channel, 
propagation in the channel improved. At 630 Hz sound was 
effectively trapped in the shallow sound channel, and propa- 
gation loss at the deep receiver was significantly greater 
than at the shallow receiver. 

Dosso and Chapman plotted relative gain between the two 
receivers according to the relationship 

Relative Gain = TL, - TL , . , (3.2) 

deep shallow 
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Figure 3.8 displays the results for a four-point range 
average centered at 24 km. The cutoff frequency of 250 Hz 
is clearly demonstrated. Below 250 Hz the SSC fails to trap 
the sound effectively, giving larger transmission losses 
than the DSC and a negative relative gain. Above 250 Hz the 
SSC traps the sound better, producing smaller transmission 
losses than the DSC and a positive relative gain. Larger 
propagation losses in the SSC below 25 Hz result from 
surface decoupling effects. 

D. MODELING PROPAGATION LOSS IN THE SSC 

Geometrical spreading is a simple model, but it is not 
able to predict the frequency and range dependence that 
Dosso and Chapman observed. A split-step parabolic equation 
(PE) model was used to calculate propagation loss based on 
the sound velocity profiles measured at ranges of zero, 16 
and 33 km along the track. [Ref. 2]. Attempts to model 
propagation using only these three profiles, however, were 
unsatisfactory. The 16 km profile indicated that the axis 
of the SSC had deepened, so that the sources were no longer 
contained in the channel. This is reflected in the increased 
propagation loss over the 5-20 km range in Figure 3.5. The 
improved propagation at and beyond 20 km indicated to Dosso 
and Chapman that the axis of the SSC must have risen to 
about 100 meters again. For this reason they used the 
profile recorded at 33 km for the 20 km range also. Dosso 
and Chapman also truncated their four sound velocity 
profiles at a depth of 600 meters, introducing an acousti- 
cally transparent bottom at that depth [Ref. 8]. This step 
had the practical effects of limiting the Fourier transform 
matrix to a manageable size and shortening computational 
time somewhat. It may also have had some effect on the 
predicted losses. 
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With these four sound velocity profiles Dosso and 
Chapman's model predicted losses with an acceptable degree 
of accuracy at 25 Hz and 630 Hz (the practical lower and 
upper limits of their model). But near the cutoff frequency 
of 250 Hz, the predicted losses were not in such good agree- 
ment with observed losses. The model did, however, predict 
strong dependence on both frequency and range as Dosso and 
Chapman observed. Figure 3.9 shows the model results 
compared to actual measurements. 

E. SUMMARY OF DOSSO AND CHAPMAN'S EXPERIMENT 

Dosso and Chapman point out that the shallow sound 
channel behaved like an acoustic waveguide with optimum 
propagation at 800 Hz, with increased losses at lower 
frequencies due to diffraction and at higher frequencies due 
to absorption and scattering [Ref. 2]. The propagation was 
found to be very sensitive to changes in the environment. 
The PE model that Dosso and C' ipman used was able to account 
for these range- dependent effects with good agreement with 
experimental results at low and near-optimum frequencies. 
At intermediate frequencies the model showed range and 
frequency dependent effects, but the predicted losses were 
greater than experimental results indicated. Dosso and 
Chapman emphasize that more accurate modeling of propagation 
losses requires a much finer sampling of the environment; 
that is, sound velocity profiles at much closer intervals 
than 16 kilometers. 
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3.1 Dosso and Chapman's sound velocity profile 

at 0 km . 
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3.2 Dosso and Chapman's sound velocity profile 

at 16 km. 
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3.3 Dosso and Chapman's sound velocity profile 
at 20 and 33 km. 
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Figure 3.4 Propagation loss vs range (Dosso and Chapman). 
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Figure 3.5 Propagation loss for 800 Hz, absorption removed 

(Dos so and Chapman) . 
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Figure 3.6 Measured loss vs frequency for three-point range 

averages (Dosso and Chapman). 
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3.7 Propagation loss' vs range for deep 
shallow receivers (Dosso and Chapman). 
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Figure 3.8 Relative gain between. deep and shallow receivers 

(Dosso and Chapman; . 
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Figure 3.9 PE model results with observations for 630 and 

25 Hz (Dosso and Chapman) . 
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IV. THE NPS PARABOLIC EQUATION MODEL 



Dosso and Chapman reported good results [Ref. 2] using 
their split-step Parabolic Equation (PE) model to predict 
propagation loss in a range- dependent environment where both 
a shallow sound channel and a deep sound channel existed. 
The first step in this study was to attempt to reproduce 
those results using the PE model available at Naval 
Postgraduate School Monterey. 

A. MODEL DESCRIPTION 

The split-step Parabolic Equation (PE) acoustic trans- 
mission loss model was chosen for study because it produces 
a numerical solution for an approximation to the wave equa- 
tion and is generally considered to provide one of the most 
accurate tools for predicting transmission loss. It can 
deal with range-dependent environments (multiple sound 
velocity profiles) and with some interaction of sound with 
the sea floor. Its drawbacks are its relatively low- 
frequency capability, inability to handle all bottom loss 
conditions, and lack of flexibility in dealing with boundary 
conditions at the sea surface (the model uses a flat 
pressure-release boundary which does not account for the 
effects of wind and waves). This PE model is implemented on 
the IBM 3033 computer at the Naval Postgraduate School (NPS) 
and is accessed through an interactive program developed and 
managed by the School's Oceanography Department. 

The FORTRAN source code for this PE model was written by 
the Acoustic Environmental Support Detachment (AESD) of the 
Naval Ocean Research and Development Activity (NORDA) during 
the 1970's. A complete description- of the physics and algo- 



37 



rithms for this model are given in [Ref. 7]. The wave equa- 
tion 

[V + K q n (r,z)] p(r,z) = 0 (4-1) 



is an elliptic equation which describes acoustic pressure as 
a function of range r and depth z. It can be approximated 
by a parabolic partial differential equation of the form 



a i 

3r 



= i 



(A + B) 



¥ 



(4.2) 



where 



A = 





and 



B 




(n 2 -l) 



(4.3) 



(4.4) 



This solution makes the assumptions that the sound is 
propagating radially outward from a source in a cylindri- 
cally symmetric medium and that the range r is very large 
compared to the size of the source and the wavelength of the 
sound: these are the cylindrical symmetry and the far-field 
assumptions. The source is further considered to be oscil- 
lating harmonically at a discrete frequency. The ocean 
surface is considered to be a flat pressure-release 
boundary, and the acoustic field vanishes in the bottom at 
the maximum depth of the Fourier transform [Ref. 7]. 

The parabolic equation accounts for diffraction and 
other full-wave effects; it also allows solution for range- 
dependent environments. The parabolic equation can be 
solved numerically using the Tappert -Hardin split-step 
Fourier algorithm to yield a solution for the entire 
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acoustic pressure field as a function of range and depth 
[Ref. 7]. 

The PE solution is best where the source and receiver 
are arranged so that the path between the two makes a small 
angle with the horizontal. At larger angles, the group and 
phase velocities of the computed solution become subject to 
errors and produce inaccurate representations in the inter- 
ference pattern. This affects the accuracy of the predicted 
transmission loss for a given range. Computing time 
increases quickly with frequency. The numerical solution is 
limited by the practical Fourier transform size, which means 
that the input environment governs the extent to which the 
model can handle higher frequencies. For practical purposes 
the PE model should be used to predict propagation loss for 
frequencies below 200 Hz along waterborne or shallow bottom- 
bounce paths [Ref. 7]. 

The PE model becomes inappropriate for situations where 
the bottom slopes steeply and where large sound speed gradi- 
ents occur [Ref. 7]. 

B. VALIDATION AGAINST DOSSO AND CHAPMAN’S DATA 

This study is concerned with waterborne paths and 
shallow angles. We are not specifically interested in 
bottom interactions. Even though the frequencies in which 
we are interested may be as high as 800 Hz [Ref. 2], we are 
willing to use computational time to examine a full-wave 
solution. Especially important for comparison with Dosso 
and Chapman's experimental data is the fact that the PE 
model can utilize a range-dependent environment. The input 
parameters described in Appendix A are derived from [Ref. 2] 
and [Ref. 8]. The resulting information and displays show 
reasonable agreement with Dosso and Chapman's work. 
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Compare Figures 4.1 through 4.5 with Figure 3.4. The 
smooth, downward trending curve represents transmission loss 
due to geometric spreading alone, calculated from equation 
3.1. The measured 1/3 octave band levels of Figure 3.4 show 
signs of smoothing because many frequencies have been added 
together, lessening the multipath interference effects that 
become apparent for discrete frequencies [Ref. 2] . 

It is apparent, however, that frequencies below the 
cutoff frequency of 250 Hz propagate poorly within the 
shallow sound channel and that higher frequencies propagate 
better. An interesting outcome here is that the 2000 Hz 
experiences less transmission loss than the experimentally 
determined optimum frequency of 800 Hz at ranges greater 
than 17.5 km. Note also that the three highest frequencies 
all show sudden steep increases in transmission loss in the 
17-20 km range and that the levels never really recover. 
This is the range at which Dosso and Chapman noted a change 
in their sound speed profiles so that their sources were no 
longer within the shallow sound channel. 

Dosso and Chapman observed an increase in transmission 
loss over roughly the same ranges (Figure 3.4), followed by 
better propagation at ranges beyond 20 km. They were able 
to model these results successfully by adding another sound 
velocity profile at 20 km, using the same profile they had 
measured at 33 km. With the same input parameters the NPS 
PE model yielded a much steeper drop and did not indicate a 
return to levels expected from the experimental data. 

Dosso and Chapman's Figure 3.5 compares to Figure 4.6. 
For Figure 3.5 Dosso and Chapman plotted measured data, with 
losses due to chemical absorption removed. This figure 
clearly shows the close correspondence between the 
measurements for the optimum frequency of 800 Hz and the 
losses predicted for geometrical spreading alone. The only 
anomalies arise over the 8-20 km range, as discussed above. 
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The NPS model allows the user to select whether attenuation 
will be considered. Figure 4.6 represents the NPS PE model 
prediction for 800 Hz without attenuation. Again the steep 
drop in transmission occurs between 15 and 20 km. A compar- 
ison with Figure 4.2 shows that the two curves are essen- 
tially similar, but attenuation accounts for a 1-2 dB 
greater loss in Figure 4.2. 

Dosso and Chapman's experimental data yielded good 
results when propagation loss was plotted as a function of 
frequency for 24 and 33 km (Figure 3.6). Recall that Dosso 
and Chapman measured 1/3 octave band levels and then calcu- 
lated three-point range averages to obtain each curve. The 
NPS PE model produces transmission loss predictions for 
discrete frequencies. To obtain transmission loss for the 
desired ranges, values were read from the computer-generated 
PE model graphic output at the appropriate ranges and at one 
kilometer on either side. The three values for each range 
were averaged to produce the values used in plotting Figure 
4.7. Again multipath interference effects can account for 
the fluctuations with range and frequency. While an optimum 
frequency of 800 Hz is easy to read from Figure 3.6, it is 
more difficult to state an optimum frequency from Figure 
4.7. Frequencies in the 500-700 Hz range seem to have 
propagated best. 

Dosso and Chapman compared results from shallow sound 
channel and deep sound channel receivers in Figure 3.7. 
Figures 4.8 through 4.10 present the corresponding results 
from the NPS PE model. As expected for a frequency well 
above cutoff, the 630 Hz curves (Figure 4.8) show better 
propagation where both source and receiver were within the 
shallow sound channel, but much greater loss (up to 10-20 
dB) where the receiver was in the deep sound channel. The 
two curves are in closer agreement for 315 Hz (Figure 4.9) 
as would be expected for a frequency near cutoff, where 
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trapping and leaking are better balanced. There is still a 
10 dB better propagation for the 15-17.5 km range when both 
source and receiver are in the channel. Figure 4.10 demon- 
strates significant transmission loss for 160 Hz, which is 
well below cutoff. In this case the deep receiver experi- 
ences better propagation by 5-10 dB over almost the entire 
range, indicating that the shallow sound channel failed to 
trap this frequency effectively. 

It should especially be noted in Figures 4.8 and 4.9 
that anomalous propagation is occurring in the 17-25 km 
range. This may be due to the change in . depth of the 
shallow sound channel that Dosso and Chapman noted and tried 
to correct by adding the sound velocity profile at 20 km. 
This compensation was apparently much less successful with 
the NPS PE model than with Dosso and Chapman's model. 

Figure 4.11 corresponds to Figure 3.8. Relative gain is 
calculated according to equation 3.2, following Dosso and 
Chapman's procedure. Their results were good for 24 km, 
showing a sudden improvement above the 250 Hz cutoff 
frequency when both source and receiver were in the SSC. 
The curve for the NPS PE model at the 24 km range does not 
produce good results, probably due to the failure of the 
added profile at 20 km to compensate for the changing depth 
of the SSC. The 15 km curve, however, is taken before the 
SSC deepened and indicates better results. The 250 Hz 
cutoff frequency is clearly indicated. The NPS results, 
however, show unexpected peaks at 40-50 Hz and at 100 Hz, 
where the shallow sound channel receiver should not be doing 
as well as the deep receiver. Other frequencies below 
cutoff behave as expected. 

Together Figures 4.12 through 4.15 correspond to Figure 
3.9, where Dosso and Chapman plotted their PE model 
predictions along with experimentally measured data. Note 
that Dosso and Chapman's 630 Hz curve shows a minimum at 
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about 18 km, followed by better propagation at longer 
ranges. Their predictions correspond rather well to the 
measurements. On the other hand, Figure 4.12 predicts 
several dB better propagation for 630 Hz than was observed 
for ranges less than 15 km and several dB worse than 
observed at ranges beyond 20 km. Again the compensation 
using the 20 km profile was less successful for the NPS PE 
model than for Dosso and Chapman's model. The steep drop 
and failure to recover in the transmission loss curves are 
even more apparent in Figure 4.13 for 400 Hz and Figure 4.14 
for 250 Hz. In Figure 4.15 for 25 Hz the observed values 
match quite well with predicted values out to about 15 km. 
At that point the prediction goes wrong again and predicted 
values fall 10-15 dB below measured values. 

C. DISCUSSION OF NPS PE RESULTS 

Discounting the multipath effects on transmission loss 
which produce extreme fluctuations in transmission loss with 
both range and frequency, it appears valid to use the NPS 
Parabolic Equation model to estimate optimum frequency for 
sound propagation in a shallow sound channel. The plot of 
transmission loss versus frequency (Figure 4.7) shows an 
optimum frequency range of 500-700 Hz. This is acceptable 
in comparison with Dosso and Chapman's measured optimum 
frequency of 800 Hz for the same set of sound velocity 
profiles. This is also good enough for active and passive 
sonar operators to use in selecting the best frequency 
ranges for their equipment . 

The NPS PE model was not, however, a good estimator for 
absolute levels of transmission loss. The NPS PE curves 
showed large differences in transmission loss levels when 
compared to Dosso and Chapman's measurements. For relative 
gain between deep and shallow receivers the NPS PE model 
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demonstrated an appropriate cutoff frequency around 250 Hz, 
but with unexpected peaks at frequencies well below cutoff. 

Some of these anomalous results from the NPS PE model 
may arise from using full-depth sound velocity profiles. 
Reasoning from ray-tracing theory, a fully absorbing bottom 
placed at 2500-2900 meters may have allowed sufficient depth 
for some deep propagation paths to be refracted back toward 
the receiver. The fully absorbing bottom at 600 meters 
which Dosso and Chapman used in their modeling may have 
eliminated propagation by these paths. Interference effects 
between these paths and the ducted propagation may account 
for some of the differences between Dosso and Chapman's 
results and the NPS PE results. 
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Figure 4.1 Transmission loss vs range, 2000 Hz 

(NPS PE model). 
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Figure 4.2 



Transmission loss vs range, 



800 Hz (NPS PE model) . 
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Figure 4.3 



Transmission loss vs range, 



250 Hz (NPS PE model) . 
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Figure 4.4 



Transmission loss vs range, 



100 Hz (NPS PE model) . 
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Figure 4.5 



Transmission loss vs range, 



25 Hz (NPS PE model) . 
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Figure 4 . 6 



Transmission loss for 800 Hz, 
(NPS PE model). 



absorption removed 
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Figure 4.7 Transmission loss vs frequency (NPS PE model). 
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Figure 4.8 



Shallow vs deep receiver, 



630 Hz (NPS PE model). 
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Figure 4.9 



Shallow vs deep receiver, 



315 Hz (NPS PE model) . 
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Figure 4 . 10 



Shallow vs deep receiver, 



160 Hz (NPS PE model) . 
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Figure 4.11 Relative gain between deep and shallow receivers 

(NPS PE model) . 
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Figure 4 . 12 



Transmission loss vs range, 
(NPS PE model). 



630 Hz 
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Figure 4.13 



Transmission loss vs range, 
(NPS PE model ) . 



400 Hz 
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Figure A . 14 



Transmission loss vs range, 
(NPS PE model). 



250 Hz 
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Figure 4.15 Transmission loss vs range, 25 Hz 

(NPS PE model) . 
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V. TWO HIGH- LATITUDE SHALLOW SOUND CHANNELS 



A. OPTIMUM FREQUENCIES IN THE SHALLOW SOUND CHANNELS 

Since the NPS PE model yielded acceptable estimates of 
the optimum frequency for Dosso and Chapman’s shallow sound 
channel, the next step in this study was to apply the model 
to some other shallow sound channels. Appendix A contains 
the NPS PE model input information for two sound velocity 
profiles taken in a deep, high- latitude ocean basin during 
summer conditions (10 and 28 June). Figures 5.1 and 5.2 
show the SVPs . Note the presence of a surface layer and the 
absence of a deep sound channel. 

The SSC for 10 June is 110 meters thick. Using equa- 
tions 2.1 and 2.2 we would expect cutoff frequencies of 148 
and 173 Hz for this channel. The SSC for 28 June is 105 
meters thick, yielding cutoff frequencies of 159 and 186 Hz. 

A plot of transmission loss versus frequency is the best 
tool for determining optimum frequency. For these two 
profiles the source was modeled at 98 meters and the 
receiver at 100 meters depth, providing for propagation 
within the shallow sound channel. Again the transmission 
loss levels at each range were calculated from the NPS PE 
model output graphics, by reading levels at the center 
ranges and one kilometer on either side, then computing a 
straight arithmetic average of the three readings. In 
Figures 5.3 and 5.4 it is much more difficult to identify an 
optimum frequency for the two shallow sound channels. 

For the 28 June case (Figure 5.3) the transmission 
losses improve gradually and show optimum frequencies in the 
600-900 Hz range, but no readily identifiable effect in the 
vicinity of the cutoff frequency. Propagation at 
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frequencies below cutoff is much better than would have been 
expected from Dosso and Chapman's results. Figures 5.5 
through 5.8 are the NPS PE model results for several 
frequencies in the optimum range. Note that there is little 
difference in the overall transmission loss levels. 

The 10 June case (Figure 5.4) is even more difficult to 
analyze. Again there is no clear indication of cutoff 
frequency in the graph. Note also that the curves for 15 
and 24 km are virtually flat for frequencies below 200 Hz. 
The two curves show that transmission loss is highly depen- 
dent upon frequency, with improved transmission at 315-400 
Hz and at several frequencies in the 600-1100 Hz range. 
Surprisingly, the transmission loss for 800 Hz on the 24 km 
curve is some 3 dB greater than for 315 Hz. Figures 5.9 and 
5.10 are the NPS PE model results for the frequencies near 
optimum, and again there is not much apparent difference in 
the overall transmission loss levels. 

In an attempt to identify the optimum frequency for the 
10 June case more clearly, the investigator then introduced 
an artificially shallow bottom in hopes that this would 
simplify the propagation paths. The 10 June sound velocity 
profile was truncated at the lower boundary of the shallow 
sound channel, with a fully absorbing bottom placed at 150 
meters. Source and receiver remained at 98 meters and 100 
meters respectively; no other input parameters were changed. 
Figure 5.11 displays the results. There is still no indica- 
tion of cutoff frequency, but two of the curves do show 
gradual improvement in transmission with increasing 
frequency. The 15 and 24 km curves show smallest transmis- 
sion losses in the 700-900 Hz range. Unfortunately the 33 
km curve shows smallest transmission loss at 600 Hz and 
large loss at 700 Hz. Figures 5.12 through 5.14 are NPS PE 
model results for frequencies near optimum. Using an artif- 
ically introduced shallow bottom appears to have improved 
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the optimum frequency prediction somewhat for the 10 June 
case . 

Putting all three cases together and looking for a 
trend, we can tentatively identify a range of optimum 
frequencies around 600-900 Hz. The results for these 
shallow sound channels are not as clear or conclusive as 
Dosso and Chapman's results, but they should be adequate for 
the fleet sonar operator to use in setting up search plans. 

B. OPTIMUM FREQUENCY FOR PROPAGATION ACROSS THE LAYER 

The second part of this investigation involved the same 
sound velocity profiles with the source still located at 98 
meters near the SSC axis, but with the receiver modeled at 
20 meters, well within the surface layer. Figures 5.15 
through 5.17 show the results. 

Figure 5.15 for 28 June shows peaks at 125 and 600 Hz in 
the transmission loss versus frequency curves for 15, 24 and 
33 km. Again the concept of cutoff frequency does not 
appear to have much meaning in these cases. Figure 5.16 
shows peaks at 160-200 Hz and at 700 Hz. For the truncated 
10 June profile, Figure 5.17 shows a peak at 75-100 Hz but 
not at higher frequencies. 

Figures 5.18 through 5.20 show relative gain for our 
three cases, calculated according to the equation 

Relative Gain = TL. - TL (5.1) 

layer ssc 

These high- latitude summer profiles do not show a cutoff 
frequency between the SSC and the surface layer. 
Propagation is always better when both source and receiver 
are contained within the shallow sound channel than when the 
source is in the SSC and the receiver is in the surface 
layer. 
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Figure 5.1 High- latitude sound velocity profile, 10 June. 
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Figure 5.2 High- latitude sound velocity profile, 28 June. 
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Figure 5.3 



Transmission loss versus frequency, 28 
(NPS PE model) . 
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Figure 5.4 



Transmission loss versus frequency, 10 
(NPS PE model) . 
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Figure 5.5 



Transmission loss versus range, 
(NPS PE model) . 



28 June, 600 Hz 
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Figure 5.6 



Transmission loss versus 
(NPS PE model). 



range, 28 June, 700 Hz 
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Figure 5.7 Transmission loss versus range, 28 June, 800 Hz 

(NPS PE model). 
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8 Transmission loss versus range, 
(NPS PE model). 



28 June, 900 Hz 
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Fieure 5.9 Transmission loss versus range, 10 June, 700 Hz 
& (NPS PE model) . 
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Figure 5.10 Transmission loss versus range, 10 June, 800 Hz 

(NPS PE model). 
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Figure 5.11 Transmission loss versus frequency, 
truncated profile (NPS PE model). 
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Figure 5.12 
truncated 



Transmission loss versus range, 
profile, 700 Hz (NPS PE model). 
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Figure 5 . 13 
truncated 



Transmission loss versus range, 
profile, 800 Hz (NPS PE model). 
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Figure 5 . 14 
truncated 



Transmission loss versus range, 
profile, 900 Hz (NPS PE model). 
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Figure 5.15 Transmission loss between surface layer and SSC, 

28 June (NPS PE model). 
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Figure 5.16 Transmission loss between surface layer and SSC, 

10 June (NPS PE model). 
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j.17 Transmission loss between surface layer and SSC, 
truncated profile (NPS PE model). 
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Figure 5.18 Relative gain between surface layer and SSC, 

28 June (NPS PE model). 
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Figure 5.19 Relative gain between surface layer and SSC, 

10 June (NPS PE model) . 
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Figure 5.20 Relative gain between surface layer and SSC, 
truncated profile (NPS. PE model). 
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VI. THE NPS FACT 9H MODEL 



A. MODEL DESCRIPTION 

The Fast Asymptotic Coherent Transmission loss model, 
version 9H (FACT 9H) , like the split-step Parabolic Equation 
model, is accessed through the Naval Postgraduate School's 
IBM 3033 computer. FACT currently provides the acoustic 
predictions for the Integrated Command ASW Prediction System 
(ICAPS) on board U. S. Navy aircraft carriers and ASW 
Operations Centers (ASWOcIs) worldwide. It is also the basis 
for a number of acoustic products routinely provided by the 
Fleet Numerical Oceanography Center, Monterey, to fleet 
users. Since the FACT model and its products are so widely 
distributed in the fleet, it is of interest to know how well 
FACT performs for the shallow ound channel case. 

The FACT 9H model as insta led at the Naval Postgraduate 
School was developed for the Naval Ocean Research and 
Development Activity (NORDA) [Ref. 9]. It uses ray tracing 
theory with special corrections to improve the treatment of 
caustics, which would produce calculations indicating exces- 
sive acoustic energy at certain ranges if not corrected. 
Unlike the PE model, FACT 9H accepts only one sound velocity 
profile, located at range zero. It makes no allowance for 
bottom topography changing over the length of a track. The 
user may indicate the absorption characteristics of the 
bottom, or specify a table of values for bottom loss versus 
grazing angle. 

In addition, a semi-empirical expression is used to 
account for propagation in a surface duct . In this module 
the user may indicate the wave height, unlike the NPS PE 
model which makes no provisions for boundary conditions at 
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FACT 9H is not 



the surface. As a ray- tracing model, 
expected to deal with waveguide phenomena such as diffrac- 
tion and leakage as successfully as the PE model does. The 
surface duct module is intended to remedy this problem some- 
what and to improve the acoustic prediction in the case of 
surface ducts. No correction has yet been applied for 
shallow sound channels. 

The investigation of the FACT 9H model proceeded along 
essentially the same lines as that of the PE model. The 
FACT 9H model was first tested against the experimental 
results of Dosso and Chapman. Then the model was applied to 
the 10 June and 28 June sound velocity profiles. Input 
parameters are described in Appendix A. 

B. VALIDATION AGAINST DOSSO AND CHAPMAN'S DATA 

Since FACT 9H does not allow for a range- dependent envi- 
ronment or for wave e: fects like diffraction and leakage, it 
was not expected to p< :form very well for the shallow sound 
channel observed by Dosso and Chapman. Again the first test 
applied to the FACT 9H output was a plot of transmission 
loss versus frequency (Figure 6.1). As for the PE model, 
the levels were obtained by taking three-point range aver- 
ages centered at 24 km. Since only the zero range sound 
velocity profile was used, there was no need to be concerned 
about the effects of the shallow sound channel axis deep- 
ening over the eight to twenty kilometer range. 

Surprisingly, the transmission loss when both source and 
receiver were within the SSC showed a gradual decrease and 
optimum values at 800-900 Hz, consistent with Dosso and 
Chapman's actual measurements. Over the entire frequency 
band tested the FACT 9H model predicted losses between 70 
and 78 dB , while Dosso and Chapman observed losses between 
about 80 and 110 dB . The transmission loss in the optimum 
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frequency range is some 10 dB less than Dosso and Chapman's 
measurements, but at least the optimum occurred at the 
proper frequencies. The transmission loss curve for the 
deep sound channel receiver showed only a fairly consistent 
increase with higher frequencies. 

Figures 6.2 and 6.3 are the FACT 9H transmission loss 
curves for the three frequencies near optimum. The shallow 
receiver curves (Figure 6.2) show more fluctuations, prob- 
ably due to multipath arrivals of sound reflecting from the 
flat ocean surface. The deep receiver curves (Figure 6.3) 
are smoothed by the partial absorption of sound in the ocean 
bottom. 

The plot of relative gain (Figure 6.4) reiterates the 
information in Figure 6.1 in a way that should indicate the 
presence of a cutoff frequency. We already know that the 
channel should exhibit a cutoff frequency at about 250 Hz, 
but there is no evidence for this in Figure 6.4. The 
shallow receiver always shows better propagation than the 
deep receiver. This result is not surprising in that the 
FACT 9H model is not designed to handle such wave effects. 

C. THE HIGH- LATITUDE SHALLOW SOUND CHANNELS 

With good predictions for optimum frequency in the Dosso 
and Chapman case, the investigation proceeded to the 10 June 
and 28 June high- lat itude , deep ocean sound velocity 
profiles. For the 28 June case (Figure 6.5) transmission 
loss in the shallow sound channel did indeed show a gradual 
rise and fall. The best frequency, however, was in the 
vicinity of 400 Hz, much lower than the 600-900 Hz predicted 
by the PE model (Chapter 5, Section A). As one might 
expect, transmission losses between the SSC and the surface 
layer are always significantly greater than when the source 
and receiver are both contained in. the SSC. 
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As for the 10 June case, FACT 9H results showed a slight 
but steady increase in transmission loss over the entire 
frequency range (Figure 6.6). The best frequency was, in 
fact, the lowest tested: 75 Hz. The transmission loss 
curve for the receiver located in the surface layer is only 
a few dB greater than the SSC curve. Recall that the NPS PE 
model also failed to produce satisfactory results for the 10 
June case. Truncating the sound velocity profile at 150 
meters with a fully absorbing bottom seemed to improve the 
PE model prediction somewhat , but Figure 6 . 7 does not 
present much of a difference for the FACT 9H model. 

Perhaps the 10 June case is a "weaker" shallow sound 
channel for modeling purposes. Dosso and Chapman's sound 
velocity profiles show a velocity gradient between 0.1436/s 
and 0.2071/s between the bottom of the surface layer and the 
axis of the shallow sound channel. The 28 June sound 
velocity profile has a velocity gradient of 0.1871/s. The 
10 June sound velocity profile, however, has a velocity 
gradient of only 0.0737/s. There may be a critical value 
for the velocity gradient between the surface layer and the 
SSC axis that will produce good results in the models. 

For now, however, it does not appear that the FACT 9H 
model is successful in predicting optimum frequencies for 
propagation of sound in shallow sound channels. In view of 
the extensive usage of the FACT model to provide acoustic 
transmission loss products to fleet operations, this problem 
warrants further attention. 
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Figure 6.1 Transmission loss versus frequency, 
Dosso and Chapman (NPS FACT 9H model). 
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Transmission loss in 
Dos so and Chapman (NPS 



the shallow sound 
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Figure 6.3 Transmission loss between SSC and DSC, 
Dosso and Chapman (NPS FACT 9H model). 
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Figure 6.4 Relative gain, Dosso and Chapman 
(NPS FACT 9H model). 
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Figure 6.5 Transmission loss versus frequency, 28 June 

(NPS FACT 9H model). 
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Figure 6.6 Transmission loss versus frequency, 10 June 

(NPS FACT 9H model). 
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Figure 6.7 Transmission loss versus frequency, 
truncated profile (NPS FACT 9H model). 
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VII. CONCLUSIONS 



A. RESULTS AND DISCUSSION 

Two acoustic transmission loss models were examined to 
determine how well they predicted the optimum frequency for 
propagation of sound in shallow sound channels. The split- 
step Parabolic Equation (PE) model performed satisfactorily, 
while the Fast Asymptotic Coherent Transmission (FACT 9H) 
model was less satisfactory. 

The PE model has two advantages: it allows for a range- 
dependent environment and it is a full-wave approximation 
for the wave equation. These two factors enable the PE 
model to deal with wave phenomena like diffraction and 
ducting with a degree of success. The PE model predicted 
optimum frequencies of 600-900 Hz for the shallow sound 
channels under study, consistent with the 800 Hz actually 
measured by Dosso and Chapman. The model seemed to produce 
better results when the sound velocity gradient between the 
overlying surface layer and the shallow sound channel was 
pronounced. The 28 June case, where the gradient was 
0.1871/s, yielded 600-900 Hz for optimum frequency, but the 
10 June case, where the gradient was only 0.0737/s, did not 
produce an optimum frequency. The 10 June prediction was 
improved somewhat by introducing a fully absorbing bottom at 
the lower boundary of the shallow sound channel. 

The FACT 9H model has its advantages, too: it is effi- 
cient in terms of computational time and it has had signifi- 
cant operational use. As a ray theory model, it was not 
expected to handle situations involving surface ducts or 
shallow sound channels as successfully as the PE model, even 
with the semi-empirical surface duct module built into this 
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version of FACT. The Dosso and Chapman sound velocity 
profiles resulted in a surprisingly good prediction of 
800-900 Hz, which agreed well with the experimental results. 
The 28 June case, however, predicted an optimum of about 400 
Hz; and the 10 June case did not predict any reasonable 
optimum. Again the strength of the sound velocity gradient 
may explain these different results. 

Of our two choices, the PE model, despite its long 
computation time, appears to be the better candidate for 
predicting optimum frequency of propagation in a shallow 
sound channel. The FACT model, despite its widespread use in 
the fleet, may not be a good predictor for optimum 
frequency. The predictions for both models seem to be 
better when there is a large sound velocity gradient between 
the surface layer and the axis of the shallow sound channel. 

B. RECOMMENDATIONS AND AREAS FOR FURTHER STUDY 

The fleet user presently has no easy access to a 
Parabolic Equation model to assess optimum frequencies for 
real-time environmental conditions. While this study 
predicts an optimum range of 600-900 Hz where both source 
and receiver are contained within 80 to 110 meter thick 
channels, it would be helpful to draw up a table for other 
source-receiver depth combinations and channels of varying 
thickness . 

Several areas merit further investigation. First, this 
study uses a fully absorbing bottom throughout. This is 
unrealistic. It would be of interest to assess bottom 
losses more accurately and determine what difference, if 
any, this makes in the prediction of optimum frequency. 
Another area for investigation is the relative strength of 
shallow sound channels. A strong channel, where the sound 
velocity gradient between the bottom of the surface layer 
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and the axis of the shallow sound channel, appears to trap 
sound more strongly and may produce a clearer prediction of 
optimum frequency. It would be helpful to find out how much 
of a velocity gradient is required to produce an optimum 
frequency from each of the models , and whether the predicted 
optimum varies with the strength of the channel. 

Finally, other acoustic transmission loss models are 
available. The implicit finite difference (IFD) PE model is 
advertised as being able to handle boundary conditions at 
the ocean bottom more successfully than the split-step PE 
model. The IFD model should be investigated. The most 
recent revision of the FACT model, version 10A, is nearly 
ready for implementation at the Naval Postgraduate School. 
Technical documents [Ref. 10] indicate that FACT 10A will, 
among other things, handl.e ducting situations and boundary 
conditions better than FACT 9H--this may affect performance 
for shallow sound channels. A third type of acoustic model 
is also available- -the RAYMODE model. This model may 
produce better results than FACT when applied to surface 
ducts and sound channels. The RAYMODE model has recently 
been designated as the Navy's new standard acoustic propaga- 
tion loss model, and its performance should be investigated 
for as wide a variety of environmental conditions as 
possible . 
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APPENDIX A 



INPUT PARAMETERS FOR THE NPS PE AND FACT 9H MODELS 

Table I describes the input options available to the 
user through the Naval Postgraduate School's split-step 
Parabolic Equation (PE) and Fast Asymptotic Coherent 
Transmission loss (FACT 9H) models. Tables II through VI 
are the sound velocity profiles used in this study. 
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TABLE I 

Input Parameters for the NPS PE and FACT 9H Models 



PARAMETERS 

Input units 

Source depths 
Receiver depths 
Frequencies 
Maximum range 
Range step 

Number of SVPs 

Depth levels 
in SVPs 

Bottom loss 



Variable 
bottom depth 

Correction for 
spherical earth 

Correction for 
for volume 
attenuation 

Vertical beam 
size of source 

Special handling 
for surface duct 

Wave height 

Critical angle 



PE MODEL 

English 
or metric 

One 

Maximum 20 
One 

Unlimited 

Minimum .01 nm 
or selected by . 
program 

Maximum 50 

Maximum 50 

1) Fully absorbing 

2) User’s loss vs 
grazing angle table 



3) User's depth vs 
attenuation table 

Maximum 100 ranges/ 
depths 

Available 

Available 

1-33 degrees 
Not required 
No 

Not required 



FACT 9H 

English 
or metric 

One 

One 

Maximum 6 

999 kilometers 

Selected by 
program 

I 

One 

Maximum 50 

1) Bottom type 
0 (reflecting) 

to 9 (absorbing) 

2) Us'r's loss 
vs gra ing angle 

t ole 

3) Internal loss 
vs grazing angle 
table 

No 

Available 

Available 

No 

Yes 

Yes 

1) Preset value 
28.65 degrees 

2) User input 
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TABLE II 

Dosso and Chapman's Sound Velocity Profile at 0 Kilometers 



DEPTH (FT) 


C (FT/ S ) 


DEPTH(M) 


C(M/S) 


0 . 


4927.82 


0 . 


1502.00 


45 . 


4927.82 


14. 


1502.00 


66 . 


4922.90 


20. 


1500.50 


87. 


4906.50 


26. 


1495.50 


115. 


4890.75 


35. 


1490.70 


137. 


4885 . 17 


42. 


1489.00 


169. 


4879.92 


52. 


1487.40 


198. 


4873.36 


60. 


1485.40 


232. 


4867.45 


71. 


1483.60 


251. 


4866.47 


76. 


1483.30 


273. 


4865.48 


83 . 


1483.00 


291. 


4864.83 


89. 


1482.80 


310. 


4864.83 


95 . 


1482.80 


325. 


4864.50 


99. 


1482.70 


338. 


4864.50 


103 . 


1482.70 


357. 


4863.52 


109 . 


1482.40 


369. 


4863.84 


112. 


1482.50 


384. 


4864.83 •• 


117. 


1482.80 


403 . 


4865.48 


123. 


1483.00 


433 . 


4866.14 


132. 


1483.20 


476. 


4867.12 


145 . 


1483.50 


510. 


4867.12 


155 . 


1483.50 


543 . 


4866.47 


166. 


1483.30 


589. 


4866.47 


179. 


1483.30 


690. 


4862.20 


210. 


1482.00 


779 . 


4859.58 


238. 


1481.20 


856. 


4858.59 


261. 


1480.90 


946. 


4857.61 


288. 


1480.60 


1026. 


4854.66 


313 . 


1479 . 70 


1115 . 


4854.00 


340. 


1479.50 


1312. 


4849.41 


400. 


1478.10 


1640. 


4849.08 


500. 


1478.00 


1969. 


4849.73 


600. 


1478.20 


2625 . 


4852.36 


800. 


1479.00 


3281. 


4856.62 


1000. 


1480.30 


3937. 


4860.89 


1200. 


1481.60 


4921. 


4870.08 


1500. 


1484.40 


6562. 


4891.07 


2000. 


1490.80 


8202. 


4916.01 


2500. 


1498.40 



See Figure 3.1. 



99 



TABLE III 

Dosso and Chapman's Sound Velocity Profile at 16 Kilometers 



DEPTH (FT) 


C (FT/ S ) 


DEPTH(M) 


C(M/S) 


0 . 
16. 
33. 
66 . 
82. 


4923.55 

4923.88 

4926.51 

4924.21 

4924.21 


0 . 
5. 
10. 
20. 
25 . 


1500.70 

1500.80 

1501.60 

1500.90 

1500.90 


98. 

131. 

164. 

197. 

230. 


4898.62 
4884.18 
4879.26 

4877.62 
4873.69 


30. 

40. 

50. 

60. 

70. 


1493.10 

1488.70 
1487.20 

1486.70 
1485.50 


262. 

295. 

328. 

361. 

394. 


4867.12 
4864.50 
4865.48 
4863.84 
4865 . 16 


80. 

90. 

100. 

110. 

120. 


1483.50 
1482.70 
1483.00 

1482.50 
1482.90 


427. 

459. 

492. 

525. 

591. 


4863.52 
4860.89 
4860.56 ■' 
4861.22 
4862.86 


130. 

140. 

150. 

160. 

180. 


1482.40 

1481.60 

1481.50 

1481.70 

1482.20 


689. 
787. 
853. 
919 . 
984. • 


4861.55 
4861.55 
4860.23 
4859 .91 
4859 .91 


210. 

240. 

260. 

280. 

300. 


1481.80 

1481.80 

1481.40 

1481.30 

1481.30 


1050. 

1115. 

1148. 

1181. 

1214. 


4858.59 

4858.27 

4857.61 

4855.64 

4853.67 


320. 

340. 

350. 

360. 

370. 


1480.90 

1480.80 

1480.60 

1480.00 

1479.40 


1247. 

1280. 

1312. 

1345. 

1411. 


4852.69 

4850.72 

4851.37 

4850.39 

4850.39 


380. 

390. 

400. 

410. 

430. 


1479.10 

1478.50 

1478.70 

1478.40 

1478.40 


1444. 
1509. 
1640. 
1969. 
2625 . 


4849.73 

4849.41 

4850.39 

4850.72 

4852.69 


440. 

460. 

500. 

600. 

800. 


1478.20 

1478.10 
1478.40 
1478.50 

1479.10 


3281. 

3937. 

4921. 

6562. 

8202. 

9514. 


4856.62 

4860.89 

4870.08 

4891.07 

4916.01 

4937.01 


1000. 

1200. 

1500. 

2000. 

2500. 

2900. 


1480.30 

1481.60 

1484.40 

1490.80 

1498.40 

1504.80 



See Figure 3.2. 
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TABLE IV 



Dosso and Chapman's Sound Velocity Profile at 20 

and 33 Kilometers 



DEPTH(FT) 


C (FT/S ) 


DEPTH(M) 


C(M/S) 


0 . 


4922.57 


0 . 


1500.40 


33 . 


4921.26 


10. 


1500.00 


6 6 . 


4896.32 


20. 


1492.40 


98 . 


4885.82 


30. 


1489.20 


131. 


4878.28 


40. 


1486.90 


164. 


4875.66 


50. ’ 


1486.10 


197. 


4875.33 


60. 


1486.00 


230. 


4873.69 


70. 


1485.50 


262 . 


4871.06 


80. 


1484.70 


279. 


4866.47 


85. 


1483.30 


295. 


4864.50 


90. 


1482.70 


312. 


4863.52 


95. 


1482.40 


328 . 


4864.17 


100. 


1482.60 


344. 


4863.84 


105. 


1482.50 


361. 


4863.84 


110. 


1482.50 


394. 


4865.16 


120. 


1482.90 


427. 


4864.83 • 


130. 


1482.80 


459. 


4865.81 


140. 


1483 . 10 


492. 


4865.48 


150. 


1483.00 


558. 


4865.81 


170. 


1483 . 10 


591. 


4864.17 


180. 


1482.60 


623. 


4863.52 


190. 


1482.40 


656. 


4860.56 


200. 


1481.50 


722. 


4860.23 


220. 


1481.40 


787. 


4857.61 


240. 


1480.60 


853. 


4856.30 


260. 


1480.20 


919. 


4853.67 


280. 


1479.40 


984. 


4850.72 


300. 


1478.50 


1083. 


4849.73 


330. 


1478.20 


1148. 


4851.05 


350. 


1478.60 


1312. 


4848 . 75 


400. 


1477.90 


1640. 


4849.08 


500. 


1478.00 


1969. 


4849.41 


600. 


1478.10 


2625. 


4851.70 


800. 


1478.80 


3281. 


4855.97 


1000. 


1480.10 


3937. 


4860.89 


1200. 


1481.60 


4921. 


4870.08 


1500. 


1484.40 


6562. 


4891.07 


2000. 


1490.80 


8202. 


4916.01 


2500. 


1498.40 


8858. 


4927.82 


2700. 


1502.00 



See Figure 3.3. 
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TABLE V 



High- latitude Sound Velocity 


Profile for 


10 June 


DEPTH (FT) 


C ( FT / S ) 


DEPTH (M) 


C(M/S) 


0 . 


4796.59 


0 . 


1462.00 


131. 


4796.59 


40. 


1462.00 


164. 


4780.18 


50. 


1457.00 


328. 


4780.18 


100. 


1457.00 


443. 


4773.62 


135. 


1455.00 


492. 


4798.23 


150. 


1462.50 


820. 


4842.52 


250. 


1476.00 


1640. 


4824.47 


500. 


1470.50 


3281. 


4847.11 


1000. 


1477.40 


6562. 


4890.42 


2000. 


1490.60 


9843. 


4943.57 


3000. 


1506.80 


13123. 


4999.34 


4000. 


1523.80 



See Figure 5.1. 



TABLE VI 



High- latitude Sound Velocity 


Profile for 


28 June 


DEPTH (FT) 


C (FT/ S ) 


DEPTH (M) 


C(M/S) 


0 . 


4806.43 


0 . 


1465.00 


115. 


4804.46 


35. 


1464.40 


164. 


4775.26 


50. 


1455.50 


262. 


4771.65 


80. 


1454.40 


328. 


4765.42 


100. 


1452.50 


459 . 


4799.54 


140. 


1462.90 


656. 


4811.68 


200. 


1466.60 


1640. 


4827.43 


500. 


1471.40 


3281. 


4847.11 


1000. 


1477.40 


6562. 


4890.42 


2000. 


1490.60 


9843. 


4943.57 


3000. 


1506.80 


13123. 


4999.34 


4000. 


1523.80 



See Figure 5.2. 
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